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Abstract 

We investigate the dynamics of the pendulum suspended on the forced Dufhng oscillator. The detailed 
bifurcation analysis in two parameter space (amplitude and frequency of excitation) which presents both 
oscillating and rotating periodic solutions of the pendulum has been performed. We identify the areas with 
low number of coexisting attractors in the parameter space as the coexistence of different attractors has a 
significant impact on the practical usage of the proposed system as a tuned mass damper. 

Keywords: Duffing oscillator, pendulum, mass tune damper, coexistence of attractors, bifurcation 
analysis. 



1. Introduction 

The dynamics of the pendulum or systems containing the pendulum is probably one of the oldest scientific 
topics. Simple parametrically excited pendulum shows extremely complex behaviour P, Miles Q shows 
that for the pendulum the route to chaos leads through symmetry breaking pitchfork bifurcation and cascades 
of period doubling bifurcations. Comprehensive analytical investigation of the pendulum with different 
forcing was presented by Bryant and Miles 0, In the forced system one can consider two main 

parameters: the amplitude and the frequency of excitation. Such a bifurcation diagram for parametrically 
forced pendulum with horizontally moving point of the suspension was presented by Bishop and Clifford for 
periodic oscillations (PO) [7] and for periodic rotations (PR) In Q bifurcation analysis was extended 
to elliptic movement of suspension point both for PO and PR. The analytical investigation of oscillating 
and rota ting motions of pendulum was done using averaging, small parameter, harmonic balance, and other 
methods[ia H Q 0, [isl . [l^ [TtI . [isl . [l9| . This analysis allows understanding of the pendulum dynamics 



in the neighbourhood of the locked periodic solutions. The interesting effect can be observed when the 
symmetry of the pendulum is broken, i.e., the imperfection term is added to the potential function [20| . In 
such an asymmetric system the symmetry breaking pitchfork bifurcation disappears and the sudden decrease 
of the amplitude of the PO at the first period doubling bifurcation can be regarded as a precursor of an 
escape or problems with system's behaviour. 

The crucial point in modelling of such systems is a good approximation of damping coefficients (viscous 
and frictional damping). There are well known methods for linear systems but as far as the pendulum 

is oscillating with large amplitude, the linear approximation does not give sufficient results. This problem 
has been solved by Xu et al. (2ll. E2|. They show an efficient method to extract the damping coefficients 
(viscous damping and dry friction) from time series of freely oscillating pendulum and consider the infiuence 
of the shaker (the source of the forcing) . 

Most previous works on the dynamics of the pendulum suspended on the forced oscillator consider the 
linear oscillators. Such a system can be considered as a modification of the classical tuned mass damper 
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H,!!!]. Early works [Hi give approximate results by the method of harmonic balance in the primary 



parametric instability zone, which allows calculation of the separate regions of stable and unstable harmonic 
solutions. Further analysis allows understanding of the dynamics around primary and secondary resonances 

[13, SIMS 

[31I I . In the recent work Ikeda proposed the usage of two pendulums mounted in the same pivot 
as a tuned mass damper [s^l ■ His experimental results show good agreement with the numerical simulations 
and his model can be considered as a good alternative to one pendulum on the pivot in the design of tuned 
mass dampers. 

A good understanding of dynamics of tuned mass damper with linear base system gives possibility to 
extend investigation to systems with non-linear base. Non-linearity in considered class of systems is usually 
introduced by changing the linear spring into non-linear one [s^, [SJI or magnetorheological damper (33 |. In 
a few papers on this topic one can find an anal ytic al study of the dynamics of Dufhng - pendulum systems 



around principal and secondary resonances |3Cl |32| . I33l l35l l36l | . The main conclusion coming from the above 
mentioned papers is that non-linear spring in the base system causes enlargement of parameters range where 
pendulum can be used as a tune mass damper. 

In this paper we consider the pendulum suspended on the forced non-linear Duffing oscillator. The 
purpose of our analysis is to study the emergence and the stability of PO and PR in two parameters space: 
the amplitude and the frequency of excitation. We identify the regions with one stable periodic solution, 
several coexisting periodic solutions, quasi-periodicity, and chaotic behaviour. 

The paper is organized as follows. In Section 2 we formulate the dimensionless equations of motion. 
The possible scenarios of pendulum's destabilization are presented in Section 3. Section 4 shows two- 
dimensional bifurcation diagrams for PO, PR, as well as one- parameter continuations for representative 
values of parameters. The influence of non-linearity of spring on damping properties of the pendulum is 
studied as well. In Section 5 we show the regions in two-dimensional parameter space where one, two, or 
several coexizted attractors can be observed. Finally, in Section 6 we summarize our results. 



2. Model of the system 

The analyzed system is shown in Fig. [1] It consists of a Duffing oscillator with a suspended pendulum. 
The Dufhng system is forced by periodical excitation and moving in a vertical direction. The position of 
mass M is given by coordinate y and the angular displacement of pendulum (position of the mass m) is 
given by angle Lp. The equations of motion can be derived using Lagrange equations of the second type. The 
kinetic energy T, potential energy V , and Rayleigh dissipation D are given respectively by the following 
formulas: ^ ^ 

T — i^iM + m)y^ — mlyip sin ip + —mPip^, (1) 

y = ifciy2 + ifc2/ + m5?(l-cos(^), (2) 

D = \ciV^ (3) 

where M is mass of the Dufhng oscillator, m is mass of the pendulum, I is length of the pendulum, fci and 
^2 are linear and non-linear parts of spring stiffness, and ci is a viscous damping coefhcient of the Duffing 
oscillator. The generalized forces are given by the following formula: 

where F [t) — Fq cosvt is a periodically varying excitation with amplitude F and frequency v, M (ip) — C2<p 
is a damping momentum with damping coefhcient C2 . The damper of pendulum is located in a pivot of the 
pendulum (not shown in Figure [I}. 

One can derive two coupled second order differential equations: 
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Figure 1: Model of system. 



(M + rn)y — mlipsmip — mlLfP' cos(p + kiy + fc2j/^ + ciy — Fq cosi^i, 



mP(p — mly sin ip + mlg sin + C2</3 = 0. 



(5) 
(6) 



In the numerical calculations we use the following values of Duffing oscillator's parameters: M — 5.0 [kg], 
ki = 162.0 [^], ^2 = 502.0 [^] , ci = 3.9 [^] and the following values of the pendulum's parameters: 
TO = 0.5 [kg], I = 0.1 [m], C2 = 0.001 [^] . We neglect static deflation of mass M. 

Introducing dimensionless time t = twi , where oj^ = jj+m natural linear frequency of Dufhng oscilla- 
tor, we reach dimensionless equations: 



where 



,2-9 „ - 



M+m 



b = 



abj sin 7 — 067^ cos 7 + a; + ax^ + dix — f cos /it, 
7 — ■^isin7 + sin 7 + ^27 = 0, 
/ = 



(7) 



■jx sm7 



(M+m)ujf 



{M+m)lwj 



di = 



{M+m)uJi 



, d2 = 



The dimensionless parameters of the system have the following values: a — 0.091, b — 3.33, a — 0.031, 
d\ = 0.132 and d2 = 0.02. Both subsystems (Duffing oscillator and the pendulum) have linear resonance 
for /i = 1.0, so around this value we expect the appearance of the complex dynamics. Amplitude / and 
frequency /x of the excitation are taken as control parameters. 

3. Destabilization of pendulum 

System (O possesses three qualitatively different regimes. The first one is the regime when the oscillations 
of the Duffing system are not large enough to destabilize the pendulum. Hence, the pendulum is at stable 
steady state 7 = 0. The second regime appears when pendulum destabilizes and start oscillating. The third 
regime is characterized by the rotating motions of the pendulum. The transition from the stable quiescence 
state 7 = to oscillations can be understood from the theoretical point of view as the destabilization of the 
invariant manifold 7 = 0. Indeed, the manifold 7 = 7 = is invariant with respect to ([7]) and the dynamics 
on the manifold is described by the single Duffing oscillator: 



Xq 



diXQ — f cos fiT, 
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with an effective mass m + M . The hnear stabihty of motions on this manifold is given by the variational 
equation: 



which has form of a linear system with respect to the variation ^7, perturbed periodically by the Duffing 
xoij). Such a parametrically perturbed system is known to possess destabilization regions (parametric 
resonance) when the frequency of the perturbation xq (t) is rationally related to the frequency of the pertur- 
bation. Since in resonances zones the Duffing oscillator is usually locked 1 : 1 with the external frequency, 
we expect the oscillation regions in the parameter space close to /i w 2 (the most prominent resonance) as 
well as ^ « 1/2, 1, 2/3, etc. 

As a result of the destabilization of pendulum, PO appear, where xq{t) is locked 1 : 1 to the external 
force and 7(t) to some other ratio depending on the resonance tongue. Since after the destabilization of 
the pendulum, the emerged periodic solution still coexists with the unstable solution (7 = 0, xq{t)), we will 
call it branching bifurcation. In Section 2] we develop a two-dimensional bifurcation diagram with respect 
to /i and /. Branching bifurcations are shown as solid lines on this bifurcation diagram in Figs. [HJa,b) 
delineating the resonance tongues. 

4. Two parameters continuation 

In this section we present two bifurcation diagrams calculated in two-parameter space: amplitude / 
versus frequency /i of excitation. We focus our attention on bifurcations of the pendulum. Duffing system, 
due to excitation, is oscillating in the whole considered range of parameters. Such plots give an overview of 
system dynamics showing the most important periodic solutions, i.e., periodic solutions with significant area 
of existence. Our calculations have been performed using software for numerical continuation AutoOTp [38]. 
As the starting points in our calculations, we use the steady state with / = 0.0 and follow it for different 
values of /i detecting the bifurcations leading to different periodic motions. Moreover, in a few cases we start 
from periodic orbits calculated by the direct integration of eq. ([7]). For integration we use the fourth order 
Runge-Kutta method. The stability of periodic solutions is given by the set of Floquet multipliers (39| . 

J^.l. Oscillatory solutions in plane (/, /i) 

In Fig. mja) we show main resonances for which the pendulum is oscillating. Different colours of 
bifurcation lines indicate the borders of different resonances tongues, i.e., areas with different locking ratio 
between the pendulum and the excitation frequency. Duffing oscillator after branching bifurcation is always 
locked 1 : 1 with excitation frequency, further bifurcations can change this ratio. The natural dimensionless 
frequency of the pendulum and the Duffing oscillator are equal to one. We find the resonances with the 
following locking rations: 1 : 1 (purple line), 1 : 2 (orange line), 2 : 1 (green line), 4 : 3 (blue line), 2 : 5 
(yellow line), and 2 : 3 (red line). These borders of the resonance tongues are of two kinds. Continuous lines 
correspond to the destabilization of the pendulum at the branching bifurcation and the appearance of PO 
with oscillating pendulum, and dashed-dotted lines to saddle-node bifurcations of PO. 

Resonance 4 : 3 appears in the saddle-node bifurcation and it is stable in a large range of parameters. 
Other resonances have qualitatively similar structure with three bifurcation curves meeting in one point: 
two solid lines and one dashed-dotted. The main resonant tongues (1:1 and 2:1) come very close to the axis 
/ = due to small friction. We illustrate the corresponding bifurcation scenarios in Fig. [UJc) using one- 
dimensional bifurcation diagrams for the case of 2 : 1 resonance. Figure [21[c) shows maximum amplitude of 
the pendulum versus / for / e (0.0, 2.5) with fixed ^ = 1.9 and /i = 2.1. Sohd and dashed lines correspond 
to stable and unstable PO respectively. The black line indicates PO calculated for ^ — 1.9. As it is easy to 
see, till / — 0.479 the pendulum is in hanging down position and only the Duffing system is oscillating. Then 
through the subcritical branching bifurcation PO (which corresponds to the oscillations of the pendulum) 
appears, but this branch of the PO is unstable (continuous line on the left side of the edge of the tongue 
(Figure [2Ia) indicate this bifurcation). For / = 0.043 the saddle-node bifurcation of PO takes place and 
this branch of PO stabilizes. For / = 1.316 one can observe the symmetry breaking pitchfork bifurcation 
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Figure 2: (colour online) Bifurcation diagrams of the PO of the pendulum for system Q in two parameter space {/, fi). (a) 
Tongues of resonances: solid lines correspond to the destabilization of the pendulum and birth of PO, dashed-dotted lines 
denote saddle-node bifurcations of PO. Different colours distinguish PO with different locking ratios between the pendulum 
and the excitation, (b) Bifurcation lines where the destabilization of PO occurs via symmetry breaking, period doubling, or 
Neimark-Saker bifurcations, (c) One-dimensional bifurcation diagrams for /i = 1.9 (black lines) and /i = 2.1 (grey lines) shows 
a route to 2 : 1 resonance in more details, (d) One-dimensional bifurcation diagram (fi = 1.275) for 2 : 1 tongue, black line 
shows PO with pendulum in hanging down position and grey one with oscillating pendulum. In (c,d) solid and dashed lines 
correspond to stable and unstable PO respectively. Abbreviations: BB - branching bifurcation, NS - Neimark-Saker bifurcation, 
PF - symmetry breaking pitchfork bifurcation, SN - saddle-node bifurcation. 
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generating two asymmetric solutions. Both branches of PO originated from pitchfork bifurcation, destabihze 
in the subcritical period doubhng bifurcations, in which the new unstable PO with doubled period is born. 
This branch of PO stabilizes in the saddle-node bifurcation. Especially interesting is the situation shown 
in the enlargement in Figure [^fc). One can see the coexistence of the 4 : 1 asymmetric PO (both Duffing 
and pendulum are quadruple) first with symmetric 2:1 PO and then after symmetry breaking pitchfork 
bifurcation of 2:1 resonance with two asymmetric 2:1 PO (after this bifurcation Duffing oscillator is locked 
2:1 with excitation - the same ratio as pendulum). Further increase of / for 4:1 asymmetric PO leads through 
the period doubling scenario to chaos. This bifurcation route shows that scenario described by Miles 0, |40| 
in the considered case of system ([T]) becomes more complicated. 

In the same plot we present scenario for /i = 2.1 (grey line) where the lower equilibrium position of the 
pendulum is destabilized by the supercritical branching bifurcation for / = 0.524 (bifurcation takes place 
at continuous line on the right side of the edge of the tongue (see Figure [^fa))) and the stability of this 
solution does not change in the considered range of the amplitude of excitation /. 

The same scenarios, with division into right (only continuous line) and left (continuous and dashed- 
dotted lines) sides of the edge of the tongues are observed for other resonances (see Figure [2l^a)). In case 
of 2:1 resonance branching bifurcation is a period doubling bifurcation while for other resonances branching 
bifurcation leads to different locking ratios. 

For 2 : 1 resonance in the range fi G (1.23, 1.317) the continuous green line is below the dashed-dotted 
green line. In this area we observe the bifurcation scenario which is shown in Fig. [^d) for /i — 1.275 in 
the range / G (0.0, 2.5) versus maximum amplitude of Duffing oscillator. The black solid line shows the 
growth of the amplitude of the Duffing oscillations in the case when the pendulum is in the lower equilibrium 
position. This PO loses its stability in the saddle-node bifurcation (the pendulum persists in the equilibrium 
position) . Then we observe the branching bifurcation of the unstable PO and the appearance of new branch 
of the PO for which the pendulum is in 2 : 1 resonance with excitation (grey line). After two saddle- node 
bifurcations the branch of the PO stabilizes for /i = 0.129 (see zoom in Fig. [^Jl). Finally, for /i = 0.1378 the 
PO loses its symmetry in the pitchfork bifurcation and through the Neimark-Saker bifurcation (/i = 0.1379) 
becomes unstable. 

In Figure [5I^b) we show main destabilization scenarios for locked PO. The resonances tongues are marked 
by grey lines in the background of the plot. The Neimark-Saker bifurcations (long curve in Figure [5)Jb)) 
destroys the 1 : 1 resonant PO, i.e., above line of this bifurcation locked 1 : 1 PO does not exist. Other 
lines are connected to stability of 2 : 1 tongue. The dashed light blue line indicates the symmetry breaking 
pitchfork bifurcation and just after it we observe the period doubling bifurcation (detailed route to chaos is 
shown in zoom in Figure [^fc)). At the end of period doubling line we detect the Neimark-Saker bifurcation, 
which merges with the symmetry breaking pitchfork bifurcation and 2 : 1 resonance curve. Other PO 
presented in Figure [DJa) are stable in the whole covering range but they are accessible only for carefully 
chosen initial conditions. 

4-.2. Rotation in plane (/, /i) 

In Figure [2I^a) we show the bifurcation diagram of rotations in two-parameter (/, fi) plane. To observe 
a rotational solutions pendulum has to undergo a global heteroclinic bifurcation. It means that this class 
of solutions could not be approached by following a steady state or PO. For continuation, one has to start 
from an integrated PR orbit. Two kinds of PR can occur in a considered system similarly to parametrically 
forced pendulum. One with a constant rotational motion in one direction are called a pure rotations (8| 
and second one with change of rotation direction are termed oscillations- rotations [41]. All lines presented 
in Fig. 3(a) are saddle-node bifurcations of PR [1, [l3- We observe the following resonant continuous PR: 
three ranges of the 1 : 1 locked PR in clockwise and counter-clockwise directions (solid grey), the 1 : 3 
in clockwise and counter-clockwise directions (solid black line). Bifurcation curves of both clockwise and 
counter-clockwise directions overlap due to the symmetry 7 — > —7. The third state is a repeated sequence 
of oscillations-rotations: one rotation in clockwise and one in counter-clockwise direction (ILIR). In Figure 
[3Ib) we show the bifurcations which destroy PR. The 1 : 1 rotational motion (solid gray line (al)) is stable in 
the large range of parameters, its stability is bounded by the Neimark-Saker bifurcation (see black line (bl)) 
in the right part of Figure[2{b)). The second 1 : 1 area is bounded by the solid gray line (a2), with increasing 
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Figure 3: Bifurcation diagram of PR for system l[7} in two parameter space {f,fJ,). Lines show bifurcations to PR with different 
locking ratios (a) and their destruction (b). 
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parameters we observe the Neimark-Saker bifurcation, which destabihzes this resonance (curve (b5)). The 
last 1 : 1 PR (hne (a3)) is bounded by the period doubhng bifurcation curve (curve (b2)). The ILli? PR 
goes through the symmetry breaking pitchfork bifurcation (curve (b4)) and the period doubhng bifurcation 
(curve (b6)) finally reaches the chaotic attractor in the period doubling cascades (the detailed description 
in Figure md)). The stability of 1 : 3 PR is bounded, on right side, by the Neimark-Saker bifurcation (curve 
(b3)), where we observe emergence of quasiperiodic motion, and further transition to chaotic attractor via 
torus breakdown with increasing of the amplitude /. On the left side this resonance disappear in period 
doubling bifurcation (line (b2)). 

4-. 3. One parameter continuation 

In this section, we present one-parameter continuations of periodic solutions (PO as well as PR) versus 
the amplitude / or frequency fi of excitation. We follow periodic solutions emerging in bifurcations presented 
in the previous sections. In Fig. Ul^a-c) we plot maximum velocity 7 while in Fig. [3]Jd) maximum angular 
position 7. PO and PR presented in Figure Hl^a) are calculated for the frequency of excitation fixed to 
IJ, = 0.92 and / is varied. When / G (0.0, 0.812) we observe the oscillation of the Duffing oscillator while 
the pendulum is in the lower equilibrium position (line (1)). At the end of this interval this state looses its 
stability in the subcritical branching bifurcation and unstable PO is created. For / = 0.162 the pendulum 
oscillations stabilize via saddle-node bifurcation and further destabilize at / = 0.274 through the Neimark- 
Saker bifurcation (curve (2)). For this value of the frequency /i there appear rotational solutions. The line 
(3) corresponds to rotational resonance 1 : 3 and the line (4) indicates ILIR PR, both solutions appear 
through the saddle-node bifurcations for / = 1.11 and for / = 1.0 respectively. Line (5) shows 9 : 9 resonance 
stable in range / 6 (0.49, 0.51) (this locked solution have small area of existence, and we do not show it in 
two-dimensional bifurcations diagram) . 

In Figure IDJb) the frequency is fixed as /it = 1.2. We observe a stable steady state of the pendulum in 
the whole range of / (line (1)). In this case, the dynamics is reduced to the motion of the forced Duffing 
oscillator. We also observe (line (2)), which corresponds to the loop in two dimensional plot (see right side of 
4 : 3 resonance in Figure |31^a)). Stable 1 : 1 rotations (line (5)) in both directions appear in the first range in 
Neimark-Saker and become unstable in period doubling bifurcation ( / G (0.39, 0.40)). In the second range 
1:1 PR stabilize by the saddle-node bifurcation and become unstable through the Neimark-Saker bifurcation 
for / G (0.743, 0.938). Last two curves are PR: 1 : 3 rotation is stable from / = 0.733 (line (3)), and ILIR 
motion is stable in range / G (0.736, 1.44) (line (4)). This last ILIR PR disappears in symmetry breaking 
pitchfork bifurcation; the detailed analysis of this PR and its bifurcations is presented in Figure |H^d). 

In Figure Hfc) we show the periodic solutions for fi — 2.25. One can observe 2 : 1 resonant PO (line 
(2)), which emerges at the branching bifurcation (/ = 1.39) from stable steady state of the pendulum (line 
(1)) and stable 1 : 1 rotation (line (3)) which starts at the saddle-node bifurcation and terminates in the 
Neimark-Saker bifurcations at the ends of the interval / G (0.147, 1.75). 

The last plot (Figure[5Jd)) shows a route to chaos starting from ILIR rotations for fixed value / — 1.25 
and variable ^. For frequency /i = 0.80, we observe the saddle-node bifurcation (on line (1)), when /z = 1.218 
the PR goes through the symmetry breaking pitchfork bifurcation and two asymmetric ILIR rotation 
appears (line (2)). Further we observe the range of existence of stable two dimensional quasiperiodic solution 
(confirmed by integration), which starts and ends in the Neimark-Saker (/ — 1.249) and the inverse Neimark- 
Saker (/ — 1.268) bifurcations (on line (2)). Finally, we observe the period doubling route to chaos. We 
show only the first solution branch with doubled period (line (3)). 

4.4- Influence of non-linearity of spring 

The characteristic of spring have a significant influence on the damping properties of the considered 
systems (sol. Is^. i.e., the energy transfer from Duffing oscillator to the pendulum during resonances. In 
Figure [S] we present these properties for two pairs of parameters: first for fi = 1.0 and / = 0.49 in 1 : 1 
resonance tongue (a) and second for /i = 2.0 and / — 1.02 in 2 : 1 resonance tongue (b). Black and gray 
colours indicate respectively the maximum position of mass M - x and the maximum angle of the pendulum - 
7 as functions of the spring non-linearity a. The continuous and dashed lines indicate respectively stable and 
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Figure 4: (colour online) Continuations of periodic solutions in one parameter / for fi = 0.92 (a), fi = 1.2 (b), and/j = 2.25 (c). 
(d) shows period doubling route to chaos starting from ILliJ rotational solution for fixed / = 1.25 and variable /i. Solid and 
dashed lines indicate the stable and unstable periodic solution respectively. Abbreviations; BB - branching bifurcation, PD - 
period doubling bifurcation, NS - Neimark-Saker bifurcation, PF - symmetry breaking pitchfork bifurcation. Other changes of 
the stability take place through the saddle-node bifurcations. 
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Figure 5: Maximum amplitude of mass (black) and maximum angle of pendulum (gray) versus non-linearity of spring a in 
range a S (—1.0, 1.0) for 1:1 resonance (a) and 2:1 resonance (b). Characteristic of spring is changing from "soft" to "hard" with 
the increasing of a. Continuous and dashed lines indicate stable and unstable periodic solutions. Abbreviation NS indicate 
Neimark-Saker bifurcation. 



unstable PO. It is easy to see that for 1:1 resonance in the case of a = Unear resonance occurs (Duffing 
osciUator is reduced to hnear oscillator). For smaller values of a the amplitude of oscillations is rapidly 
decreasing for both pendulum and Duffing oscillator up to a = —0.97 where the motion of the pendulum 
stops. Similar decreasing of the amplitude of oscillations appears for positive a, but motion terminates at 
a — 4.05. For /i = 0.05 we observe the Neimark-Saker bifurcation followed by the inverse Neimark-Saker 
bifurcation for n — 0.075. Between the bifurcation we observe a stable quasiperiodic motion. In Figure 
[51^b) one observes a stable 2 : 1 PO in the whole range of a. When the system is changing from soft to 
hard characteristic of spring amplitudes of the Duffing oscillator and the pendulum are increasing. We do 
not observe a resonance of Duffing system. In 1:1 resonance we observe a decrease of the amplitude of the 
oscillation of mass M, while in others resonances (similarly as for 2:1 PO) this amplitude increases due to 
energy transmitted from the oscillating pendulum. 



5. Coexistence of solutions 

As follows from the previous sections, see e.g. Fig. Hfa-d), multiple stable periodic solutions often coexist 
for the same parameter values. In this section we illustrate basins of attraction of different stable periodic 
solution. Numerical package Dynamics 2 [i^ is used. We choose six representative sets of parameters (/, /i) 
and show the calculated basins in plane (7,7) (angular displacement and velocity of pendulum), where 
7 is mod27r. Initial conditions of the pendulum are taken in ranges: 7 G (— tt, tt) and 7 G (4, —4), initial 
conditions of Duffing oscillator are fixed for each plot and have the following values: = 0, ig = (a,b,d,e), 
xq = 2.0, xq = 5.0 (c, f). Since different initial conditions of the Duffing oscillator may lead to different 
attractors in plane (7,7), the obtained figures show two-dimensional cross-sections of the four directional 
phase space (plus the phase of the perturbation). As one can see in Fig. [3] and Fig. [7] most of the periodic 
solutions are accumulated in range fi G (0.7, 1.3), so we calculated four out of six the basins in this area. In 
Fig. [6la) (/i = 0.9 and / = 0.5) we find five attractors: steady state, pair of period nine motion and chaotic 
motion which bifurcates from 1 : 1 PO. Then in Fig. [HJb) {fj, = 1.03 and / — 1.25) one can observe five 
attractors. Two of them are symmetric pairs of rotations (1:3 resonance), two corresponds to symmetric 
1 : 1 PO and the last one is the ILIR PR. Next plot {fi ^ 1.2 and / — 0.75) includes seven attractors: two 
symmetric pairs of PR (1:1 and 1 : 3), 4 : 3 PO, equilibrium of pendulum and ILIR PR. In Fig. [DJd) 
(/i = 1.229 and / = 1.0) we detect six solutions, two of them are quasiperiodic bifurcated from symmetric 
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pair of 1 : 1 PR , pair of 1 : 3 PR, steady state of pendulum and lili? PR. For larger values of excitation we 
also find chaotic attractor, e.g. for /z = 1.6, / = 1.95 (see Fig. EJe)). Usually the chaotic solution dominates 
the whole phase space and the coexisting attractors have small basins of attraction, and such a situation is 
also observed in the investigated system ([7]). In the last Fig. we show a case where we do not observe 
fractal basins of attraction. Most of the phase space is dominated by period two symmetric PO, only in 
small range one can see basins of symmetric pair of 1:1 PR. 

Generally, in area where lili? PR exists its basin of attraction dominate in phase space. The 4 : 3 
resonance could be observed for all pairs of parameters used in Fig. Ela-d). Nevertheless, we see this 
attractor only in Fig. EIc), where different initial conditions of Duffing are used. This is the evidence that 
not only sensitivity on initial state is observed for pendulum but also for Duffing. Varying initial condition 
of Duffing moves cross-section of the phase space and changes the set of accessible attractors. 

Almost all basins of attraction have a fractal structure, so reaching the given solution is strongly de- 
pendent on initial conditions. From the practical point of view it is important to know the area with a 
small number or even one solution [43*1. In such ranges one can be sure that the system approaches the 
expected solution. In Fig. [7] we marked the areas with different types of attractors: black colour indicates 
one attractor (four locked PO, excluding 4 : 3 resonance, between two branching bifurcation lines on the 
left and right side - the edges of the resonance tongues) , grey colour refers to two coexisting solutions (the 
same as for black but with coexisting steady state of the pendulum) . In the hatched area we observe the 
coexistence of PR (two symmetric pairs of 1 : 1 or 1 : 3 or ILIR PR) and the steady state of the pendulum. 
The largest area with one attractor is a tongue of 2 : 1 resonance. For other resonances (1 : 1, 2 : 3, 2 : 5, 
1:2) these areas are small. Especially, it is surprising for 1 : 1 where the resonance tongue is large in the 
parameter space but only near its edge we do not find the second attractor. We do not mark the areas where 
only Duffing system is oscillating and pendulum is stable equilibrium position. In this case the dynamics of 
the system is reduced to the oscillations of mass (M -|- m). 

6. Conclusions 

We present a comprehensive numerical analysis of the forced Duffing oscillator with the suspended pen- 
dulum. We show two dimensional bifurcation diagrams with the most representative periodic solutions and 
demonstrate the bifurcation route to the locked resonances. The linear resonance of both subsystems is 
observed for /i = 1.0 and around this value we find complex dynamics with many coexisting attractors, 
not only periodic but also quasiperiodic and chaotic ones. In the principal resonance zone the pendulum 
oscillations decrease the oscillation amplitude of the Duffing oscillator so we observe an energy transfer from 
the Duffing oscillator to the pendulum. These damping properties are observed only for 1 : 1 oscillatory 
resonance, for other locked solutions (oscillations or rotations) the pendulum oscillations increase the os- 
cillations amplitude of the Duffing oscillator (the control and the possible decrease of the amplitude for 
other resonances will be considered in future work). Contrary to complex dynamics around 1 : 1 principal 
resonance in the neighbourhood of 2 : 1 parametric resonance, we find two large ranges in parameters space 
with only one attractor (2 : 1 locked oscillations) and symmetric pair of 1 : 1 rotations respectively. From 
the practical point of view (certainty of reaching the desired attractor) such a situation is very useful and 
rare in non-linear system. We also compare the infiuence of non-linearity of spring on the amplitude of os- 
cillations For principal 1 : 1 resonance, strong non-linearity (hardening or softening) causes lower amplitude 
of oscillations. For 2 : 1 parametric resonance only the decrease of a into the softening direction causes 
the decrease of the oscillations amplitude. The existence of fractal basins of attraction is not surprising for 
non-linear systems with the attached pendulum. Hence, we can consider only the probability of reaching 
the chosen attractor and never have a certainty where the dynamics of systems evolves. This is crucial when 
the pendulum is working as a tune mass damper so we show the areas with a low number of coexisting 
solutions in the parameters space. 
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Figure 6: (colour) Basins of attraction in the plane of the pendulum variables (7, 7) for = 0.9, / = 0.5 (a); = 1.03, / = 1.25 
(b); n = 1.2, / = 0.75 (c); ^l = 1.229, / = 1.0 (d); /i = 1.6, / = 1.95 (e), and ^l = 2.37, / = 2.25 (f). The colours for attractors 
their basins are shown in legend. 
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Figure 7: Two-parameter bifurcations diagram in plane (/, /x) with PO and PR. Black colour indicates one attractor (for locked 
PO, excluding 4 : 3 resonance, between period doubling bifurcation lines on the left and right side - edges of the resonance 
tongues), gray colour shows two coexisting attractors (the same as for black but with coexisting stable steady state of the 
pendulum). In the hatched area we observe the coexistence of stable rotations (three areas with: pairs of symmetric 1 : 1 or 
1 : 3 or ILliJ PR) and stable steady state of the pendulum. 
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